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Abstract 

Signal transducers and activators of transcription (STATs), tlieir inliibitors and cyclooxygenase-2 (COX-2) participate in 
transformations of many various types of cancers. The aim of the present study was to evaluate the relationship between 
STAT5A/B, COX-2, and PIAS3 mRNA expression and tumor staging, metastasis status, and histopathological subtype in 71 
patients with confirmed non-small cell lung cancer (NSCLC) diagnosis. Total RNA was isolated from NSCLC tissue samples 
and the expression of the studied genes was assessed using TaqMan probes in real-time PCR assay. The expression levels of 
STAT5A, STAT5B, and COX-2 genes were increased in 69%, 79%, and 71% NSCLC samples respectively, while PIAS3 expression 
was decreased in the majority (69%) of the studied tissues. Statistically significant differences were observed between STAT5 
isoforms (P = 0.0008), with higher expression of STAT5B. We found statistically significant positive correlation between 
STAT5B and COX-2 (rho = 0.045), and significant negative correlation between STAT5B and PIAS3 (rho = -0.049). The negative 
correlation between STAT5B and PIAS3 (rho= -0.43) was also observed in T2a-HT2b tumor group. Additionally, STAT5B and 
COX-2 expression levels were significantly different between Tla-nTlb and T2a-HT2b tumors (P = 0.002 and P = 0.041, 
respectively), with higher expression of both genes in T2 tumor stage. PIAS3 expression was significantly lower in NSCC 
subtype as compared with SCC subtype (P = 0.017). Also, STAT5A and STAT5B immunoexpression was assessed, and the 
results indicated significantly higher protein levels in NSCLC patients as compared with controls (P = 0.048 and P = 0.034, 
respectively). High STAT5B immunoexpression was positively correlated with STAT5B gene expression in tumors 
(rho = 0.755). STAT5B protein level was also significantly higher in T2a-HT2b tumors, reflecting high STAT5B gene expression 
in this group. There was no statistically significant association between mRNA and protein expression levels of the studied 
genes and patients' characteristics: age, gender, smoking. The obtained results highlight the importance of the genes 
STAT5B and COX-2 in lung cancer progression. 
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Introduction 

It is documented that despite of more and more modern 
treatment options, lung cancer is one of tlie leading reason of 
cancer related deatlis in the world. Non-small cell lung cancer 
(NSCLC) is recognized as the most common - accounting for 7,5- 
85% - among all lung malignant tumors [1]. With advances in 
molecular biology, the detection and analysis of changes in 
expression levels of many important genes involved in signaling 
pathways may supply predictive molecular markers harboring 
diagnostic and/or prognostic value for NSCLC. 

In many human cancers, including NSCLC, one of the key 
pathways that promote cellular survival or cell growth is Janus 
kinase /signal transducers and activators of transcription (JAK/ 
STAT) pathway. It is one of the pleiotropic cascades of molecules 
involved in signal transduction for proliferation, development and 



apoptosis [2,3]. STAT (signal transducers and activators of 
transcription) protein family of transcription factors consists of 
seven members: 1-4, 5A, 5B, and 6 [2]. Among them, and beside 
STATS, the oncogenic activity of STAT5 was documented both 
in vitro and in vivo [4,5]. After phosphorylation, the two STAT5 
proteins, STAT5A and STAT5B, as homo- or heterodimer, may 
translocate to the nucleus and bind to specific STAT5 response 
elements of target gene promoters, thus influencing their 
expression [6]. It is recognized that STAT3A and STAT5B 
isoforms are encoded by two tandemly linked genes on chromo- 
some 17qll.2 [7]. They act as independent transcription factors 
[8] and modulate important cellular processes in different manner 
in normal and malignant cells [9]. Generally, active STATS 
promotes cell cycle progression, proliferation, invasion, angiogen- 
esis, and inhibits apoptosis. 
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Table 1. Histopathological verifications of NSCLC tissue samples. 



Histopathologlcal type of NSCLC 


squamous cell carcinoma (SCC) 


41 (57.75%) 




non-squamous cell carcinoma (NSCC) 


30 (42.25%) 


AJCC* 


AJCC lA 


14 (19.72%) 




AJCC IB 


11 (15.49%) 




AJCC HA 


13 (18.31%) 




AJCC IIB 


10 (14.08%) 




AJCC IIIA/IIIB 


23 (32.39%) 


pTNM* 


TIa+TIb 


19 (26.76%) 




T2a+T2b 


33 (46.48%) 




T3+T4 


19 (26.76%) 



*AJCC - American Joint Committee on Cancer Staging according to the lASCLC Staging Project 7th ed. (2010) Cancer. 
*pTNM — post-operative Tumor Node Metastasis classification according to the WHO Histological Typing of Lung Tumor. 
doi:l 0.1 371 /journal.pone.Ol 04265.t001 



The overexpression of STAT5 has been recognized in several 
types of human tumors, mainly in breast and prostate cancers [10- 
13]. However, as so far, the data on the role of STATS in NSCLC 
cells, as well as on its activation status in NSCLC is still very 
limited. Recendy, it has been documented that besides several 
cytokines, hormones and growth factors, EOF influences the 
STATS expression in human lung adenocarcinoma cell line thus 
leading to the increased cyclooxygenase-2 (COX-2) expression 
[14]. COX-2, belonging to the COX enzyme family (COX-1, 
COX-2 and COX- 3) [15], is a key enzyme in the biosynthesis of 
prostaglandins (PG). COX-2 is involved in the initiation and 
progress of tumors in situ and its overexpression is frequently 
recognized in many tumor types, including NSCLC [16-18]. 
However, carcinogenic effect of COX-2 upregulation in relation 
to the expression level of STATS in NSCLC patients has not been 
investigated yet. 

On the other hand, as a modulator of activity of STATS, a 
protein inhibitor of activated STAT3 (PL\S3), which regulates 
dilferent DNA binding transcription factors implicated in the 
immune response (e.g., NFkB, SMAD, and MITF), was recog- 
nized in breast cancer [19]. To date, only a small number of 
reports focused on PL\S3 participation in cancers involving lung 
tumors, has been published [20-23]. Additionally, although it is 
known that PL\S3 plays a vital role in oncogenic process 
influencing STAT3 protein, interaction between PL\S3 and 
STAT5 has not been fuUy recognized yet. The aim of our study 



was to determine the relationship between STATS, PIAS3 and 
also COX-2 and their reciprocal relationship on transcriptional 
level in NSCLC patients. To achieve this goal, we assessed the 
mRNA expression of these genes and their association with 
histopathological features of NSCLC tumors as well as clinical 
characteristic of patients. The prespecified hypothesis tested was 
that STATS, COX-2 and PIAS3 expression levels were modified 
in non-small cell lung cancer, playing a role in lung carcinogenesis. 
Additionally, we analyzed the levels of STAT5A and STAT5B 
proteins in the studied samples. 

Materials and Methods 

The study has been approved by the Ethical Committee of the 
Medical University of Lodz, Poland no. RNN/64/1 1/KE. 
Written informed consent was obtained from each patient. 

1. Characterization of the NSCLC tissue samples and 
patients clinical characteristics 

Biological material (lung tissue) was obtained from 84 patients 
admitted to the Department of Thoracic Surgery, General and 
Oncologic Surgery, Medical University of Lodz, Poland, between 
July 2010^une 2012. Based on the results of preoperative 
cytological/histological assessment, the patients were qualified 
for surgery and were treated by either lobectomy or pneumonec- 
tomy. 



Table 2. Characteristics of NSCLC patients in terms of tobacco addiction and consumption (the amount of cigarettes smoked per 
day and FYs). 





Tobacco addiction and consumption 




n = 71 


Smokers 




66 (92.96%) 


Non-smol(ers 




5 (7.04%) 


Tlie smoking period 


<40 years 


37 (52.11%) 




>40 years 


29 (40.85%) 


Tlie amount of cigarettes smoked 


10-15 cigarettes per day 


6 (9.09%) 




20 cigarettes per day (1 pack) 


43 (65.15%) 




30-40 cigarettes per day (1.5-2 packs) 


1 7 (25.76%) 


Pack Years (PYs) 


up to 40 PYs 


30 (45.45%) 




=:40 PYs 


36 (54.54%) 



doi:1 0.1 371 /journal.pone.Ol 04265.t002 
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Table 3. The expression levels (RQ values) of the studied genes in individual histopathological NSCLC subtypes. 



Gene 


Histopathological NSCLC subtype 


Mean RQ value (range) 


Number (percentage) of samples with 








RQ value^l 


RQ value^l 


STATSA 


sec (n = 41) 


3.25 (0.25-4.49) 


31 (75.60%) 


10 (24.39%) 




NSCC (n = 30) 


2.155 (0.05-4.54) 


18 (60.00%) 


12 (40.00%) 




Total (n = 71) 


2.75 (0.05-4.54) 


49 (69.01%) 


22 (30.99%) 


STAT5B 


sec (n = 41) 


8.92 (0.67-59.29) 


/OC ^"70/ \ 
3D \GD.D/ /O) 


O [ 1 4.D370J 




NSCC (n = 30) 


9.26 (0.77-864.17) 


21 (70.00%) 


9 (30.00%) 




Total (n = 71) 


9.13 (0.77-864.17) 


56 (78.87%) 


15 (21.13%) 


COX-2 


sec (n = 41) 


9.34 (0.10-58.21) 


29 (70.73%) 


12 (29.27%) 




NSCC (n = 30) 


4.89 (0.03-97.53) 


21 (70.00%) 


9 (30.00%) 




Total (n = 71) 


3.17 (0.24-20.89) 


50 (70.42%) 


21 (29.58%) 


PIAS3 


sec (n = 41) 


0.21 (0.16-1.74) 


14 (34.15%) 


27 (65.85%) 




NSCC (n = 30) 


0.94 (0.02-1.12) 


8 (26.67%) 


22 (73.33%) 




Total (n = 71) 


0.54 (0.02-14.91) 


22 (30.99%) 


49 (69.02%) 



doi:l 0.1 371 /journal.pone.Ol 04265.1003 



Immediately after resection, lung tissue samples (100-150 mg) 
and the adjacent non-cancerous macroscopicaUy unchanged 
tissues (100 mg; 10 cm distant from the primary lesion) obtained 
from the same patients were placed in a stabilization buffer 
RNAlater®. Each tissue sample was divided into smaller parts (SO- 
SO mg) for individual analysis. AU samples were frozen at — 80°C. 

The resected lesions were post-operatively histhopathologically 
evaluated and classified according to the AJCC staging as well as 
TNM classification (pTNM). Based on the results, from the group 
of 84 patients, 1 2 of them were excluded due to a concomitant 
malignancy and the suspected metastasis to the lung, in the case of 
one patient, there was no full clinical information. The final study 
group included seventy one (n=71) patients with confirmed 
NSCLC diagnoses. Histopathological assessments of tumor 
specimens were obtained from pathomorphological reports, and 
were as follows: squamous cell carcinoma (SCC, n = 41), non- 
squamous cell carcinoma (NSCC, n = 30). Histopathological 
verifications of non-smaU cell lung carcinoma tissues are shown 
in Table 1. The studied group consisted of 25 women, mean age 
63±8.717 and 46 men, mean age 65±8.234. AU cases were 
primary tumors without chemo- or radiotherapy treatment. The 
smoking history was available for all patients: 5 patients were non- 
smokers, and 66 were smokers or former smokers. They were 
divided into groups according to their smoking habits: time of 
tobacco addiction and amount of cigarettes smoked - the latter 
was presented as Pack Years (PYs) and was calculated according to 
the NCI Dictionary of Cancer Terms (1 Pack Year is equal to 20 
cigarettes smoked per day for 1 year), see Table 2. 

2. RNA extraction, real-time PCR 

Total RNA was extracted from lung tissues (cancer tissue 
obtained from the center of lung lesion and macroscopicaUy 
unchanged lung tissue obtained from the most distant site from the 
resected lesion) using Universal RNA Purification Kit (Eurix, 
Poland), and according to the manufacturer's recommendations. 
The quality and quantity assessments of RNA samples were 
determined by minielectrophoreses in polyacrylamide gel (AgUent 
2100 Bioanalyzer, Agilent, USA), using RNA 6000 Pico/Nano 
LabChip kit (Agilent Technologies, USA). Complementary DNA 
(cDNA) was transcribed from 100 ng of total RNA, using a High- 
Capacity cDNA Reverse Transcription Kit (Applied Biosystems, 



USA) in a total volume of 20 |J.l per reaction. Reverse transcription 
(RT) master mix contained: lOx RT buffer, 25 x dNTP Mix 
(100 mM), lOx RT Random Primers, MultiScribe™ Reverse 
Transcriptase, RNase Inhibitor and nuclease-free water. RT 
reaction was performed in a Personal Thermocycler (Eppendorf, 
Germany) in the foUowing conditions: 10 minutes at 25°C, 
foUowed by 120 minutes at 37°C, then the samples were heated 
to 85°C for 5 seconds, and hold at 4°C. The relative expression 
was assessed using TaqMan probes: Hs00234181_ml, 
Hs00273500_ml, Hs00153133_ml, H,s00180666_ml for the 
studied genes STATSA, STAT5B, COX-2 and PIAS3, respec- 
tively, as weU as for li-actin (ACTB, Hs99999903_ml), as the 
reference gene. The PCR mixture contained: cDNA (1 to 100 ng), 
20 X TaqMan® Gene Expression Assay, 2x TaqMan® Gene 
Expression Master Mrs, RNase-free water in a total volume of 
20 [ll. The qPCR reactions were performed in Applied Biosystems 
7900HT Fast Real-Time PCR System for 39 cycles, with 
annealing temperature of 60°C, repeated 3 times for each sample. 
The relative expression of the studied samples were assessed using 
the Comparative delta-delta Cx method (TaqMan Relative 
Quantification Assay software) and presented as RQ^ value. 



E 25 
O 2 



□ Mean 
I MeanJSD 



Figure 1. Box-and-whisker plots, representing 5771 755 expres- 
sion levels (mean RQ values) in the studied tumor groups 
classified according to pTNM staging in relation to their size. 

doi:1 0.1 371/joumal.pone.01 04265.g001 
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Table 4. The results of Neuman-Keuls' multiple comparison test regarding differences in STAT5B RQ values in the studied tumor 
groups classified according to pTNM staging in relation to the tumor size. 









TIa+TIb 


T2a+T2b 


T3+T4 


T1a+T1b 




0.001782 


1000000 


T2a+T2b 


0.001782 




0.723456 


T3+T4 


1 .000000 


0.723456 







doi:l 0.1 371 /journal.pone.Ol 04265.t004 

adjusted to ji-actin expression level. Macroscopically unchanged 
lung tissue served as calibrator sample. 

3. Tissue homogenization, ELISA assays 

Lung tissue samples (30-50 mg) were rinsed in ice-cold PBS 
bufier (0.01 mol/1, pH 7.0-7.2) and homogenized in 5 ml of PBS. 
The resulting suspension was subjected to two cycles of freezing 
and thawing. Then, the homogenates were centrifuged for 
5 minutes at 5000 xg, the supernatant was removed and the 
suspension was ahquoted and stored at — 80°C until further 
analysis. 

STAT5A and STAT5B immunoexpression levels in lung tissue 
homogenates were assessed using ELISA Kit for Signal Trans- 
ducer And Activator Of Transcription 5A (SEB738Hu, Uscn Life 
Science Inc., China) and ELISA Kit for Signal Transducer And 
Activator Of Transcription 5B (SEB727Hu, Uscn Life Science 
Inc., China), according to the manufacturer's procedure. The 
intensity of the final colorimetric reaction, in proportion to the 
amount of protein bound, was measured in a plate reader 
(ELx800, BioTek) at 450 nm. The obtained results were compared 
to the standard solution of known concentrations (0.312-20 ng/ 
ml). 

4. Statistical analysis 

ANOVA Kruskal-WaUis' test and U Mann-Whitney test, were 
used to compare the levels of relative expression values (RQs) of 
STAT5A, STAT5B, COX-2 and PIAS3 among NSCLC subtypes 
(SCC, NSCC). 



>< 

8 2.5 



3 Mean 
. □ MemSEM 
TIa+TIb iMeaniSD 



Figure 2. Box-and-whisker plots, representing tKie expression 
levels (mean RQ values) of COX-2 gene in the studied tumor 
groups classified according to pTNM staging in relation to 
their size. 

doi:1 0.1 371/journal.pone.01 04265.g002 



The Spearman's rank correlation coefficient, U Mann-Whit- 
ney's test and ANOVA Kruskal-Wallis' test were performed in 
order to evaluate the relationship between the expression levels 
(Revalues) of the studied genes {STAT5A, STAT5B, COX-2 and 
PIAS3), between RQs and the examined parameters (patient 
characteristics: age, gender and tumor staging according to pTNM 
and AJCC classification), as well as to find correlations regarding 
STAT5A and STAT5B protein levels in the studied samples. The 
accepted level of statistical significance was estimated at P<0.05. 
The results of relative expression levels of the studied genes are 
presented as RQ, (relative quantification) means ± SEM and RQ 
means ± SD values. Statistica for Windows 10.0 program was 
apphed for calculations. 

Results 

1 . Relative expression levels of STAT5A, STAT5B, COX-2 
and PIAS3 

Relative expression levels of the studied genes {STAT5A, 
STAT5B, COX-2, PIAS3), expressed as RQ, values, were 
determined using delta-delta Cx method, adjusted to the 
expression of ji-actin (endogenous control) and relative to the 
expression level of calibrator (macroscopically unchanged lung 
tissue) for which RQ,= 1. 

The obtained RQ_ values for the individual genes were 
correlated with histopathological NSCLC subtypes (SCC, NSCC), 
tumor staging (pTNM, AJCC), patients' age, gender and smoking 
history. 

STAT5A expression was increased (RQ, value>l) in the 
majority (69%) of the studied samples. Regarding the individual 
histotypes, it was increased in a greater proportion in SCC group, 
but there were no statistically significant differences. STAT5B 
expression was increased in all studied histopathological NSCLC 
subtypes, in a range of 70%-85%, depending on a histotype. 
Expression levels of STATS isoforms {STAT5A and STAT5B) 
were statistically significant different (P = 0.0008; Spearman's rank 
correlation coefficient), with higher expression level displayed by 
STAT5B. The expression level of COX-2 was nearly uniformly 
increased in all NSCLC subtypes, at the level of 70%. In contrast, 
PIAS3 expression was decreased (RQ, value<l) and it was 
observed in all studied histopathological NSCLC subtypes, in a 
range of 66%-73%, depending on a histotype. Statistically 
significant differences were found between SCC and NSCC 
subtypes (P = 0.017; U Mann- Whitney's test), with significandy 
lower gene expression observed in NSCC samples. The results 
regarding genes' expression levels are summarized in Table 3. 

2. Statistical analysis of relationship between studied 
gene expression levels and clinical features of patients 
and tunnor characteristics 

There was no statistically significant correlation between RQ, 
values of STAT5A and the clinical features of NSCLC patients. 
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Table 5. The results of Neuman-Keuls' multiple comparison test regarding differences 
groups classified according to pTNM staging in relation to the tumor size. 


in COX-2 RQ values in the studied tumor 




Tia+tib T2a+T2b 


T3+T4 


T1a+T1b 0.02459 


1 000000 


T2a+T2b 0.02459 


0.27776 


T3+T4 1 .000000 0.27776 



doi:l 0.1 371 /journal.pone.Ol 04265.t005 



i.e., patients' age (P = 0.48), gender (P=0.54), and liLstory of 
smoking assessed as PY (P = 0.87), as well as histopathological 
features of tumor, i.e., pTNM classification (P = 0.34), AJCC 
classification (P = 0.83), NSCLC subtypes (SCC vs NSCC, 
P = 0.62) (ANOVA Kruskal-Wallis test, U Mann-Whitney's test 
followed by Spearman's rank correlation coefficient). 

Similarly, statistical analysis concerning STAT5B did not reveal 
significant correlations between its expression levels and patients' 
age (P = 0.59), gender (P = 0.45), history of smoking (P = 0.66), 
NCSLC subtypes (P = 0.85) and tumor AJCC staging (P = 0.23) 
(ANOVA Kruskal-WaUis test, U Mann- Whitney's test followed by 
Spearman's rank correlation coefficient). However, statistically 
significant differences were found in relation to tumor size 
classified according to pTNM staging (P = 0.034; ANOVA 
Kruskal-Wallis test), as presented in Figure 1. Neuman-Keuls' 
multiple comparison test revealed statistically significant differenc- 
es between Tla+Tlb vs T2a+T2b group (see Table 4). 

There was no statistically significant correlation between RQ^ 
values of COX-2 and patients' age (P = 0.09), gender (P = 0.44) and 
history of smoking (P = 0.34) (ANOVA Kruskal-Wallis test, U 
Mann- Whitney's test followed by Spearman's rank correlation 
coefficient). Similarly, no statistically significant correlations were 
found between RQ^ values of COX-2 and histopthological NSCLC 
subtypes: SCC vs NSCC (P = 0.5 1) and tumor AJCC classification 
(P = 0.28) (U Mann-Whitney's test). However, the association 
between tumor size according to pTNM staging and COX-2 
expression level was statistically significant (P = 0.041; ANOVA 
Kruskal-Wallis test), as presented in Figure 2. Neuman-Keuls' 




Figure 3. Box-and-whisker plots, representing /V>455 expression 
levels (mean RQ values) in the studied tumor groups classified 
according to the histopathological NSCLC subtypes. 

doi:1 0.1 371 /journal.pone.Ol 04265.g003 



multiple comparison test revealed statistically significant differenc- 
es in RQ^ values between Tla+Tlb and T2a+T2b group (see 
Table 5). 

Regarding PI AS 3 expression, statistical analysis did not reveal 
any significant correlations between the RQ_ values and patients' 
age (P = 0.48), gender (P = 0.37), and history of smoking (P = 0.58) 
(ANOVA Kruskal-WaUis test, U Mann- Whitney's test followed by 
Spearman's rank correlation coefficient), as well as tumor staging 
according to pTNM (P = 0.22) and AJCC classification (P = 0.12) 
(ANOVA Kruskal-Wallis test). Statistically significant differences 
were found between histopathological NSCLC subtypes 
(P = 0.017; U Mann- Whitney's test), as presented in Figure 3. 

3. Statistical analysis of the reciprocal relationship 
between the expression of the studied genes 

Finally, we assessed the reciprocal relationship between the 
expression levels of the studied genes. Spearman's rank correlation 
coefficient revealed statistically significant negative correlation 
between STAT5B and PIAS3 genes (rho= -0.049, P = 0.04; 
Spearman's rank correlation) and positive correlation between 
STAT5B and COX-2 genes (rho = 0.045, P = 0.027; Spearman's 
rank correlation). Additionally, statistically significant negative 
correlation between STAT5B and PIAS3 genes was observed in 
T2a+T2b tumor group (rho = —0.43, P = 0.041; Spearman's rank 
correlation). 

4. STAT5A and STAT5B innmunoexpression levels 

A quantitative measurements of STAT5A and STAT5B in 
NSCLC and control samples were performed using ELISA 
method. The mean values (ng/ ml) obtained in the studied samples 
are presented in Table 6. Statistical analysis revealed that both 
STAT5A and STAT5B protein immuiioexpression levels were 
significandy higher in lung cancer samples as compared with 
normal (macroscopicaUy unchanged, control) lung tissues 
(P = 0.048 in case of STAT5A and P = 0.034 in case of STAT5B; 
U Mann-Whitney's test). However, the differences between 
NSCLC histopathological subtypes (SCC vs NSCC) weren't 
statistically significant, neither for STAT5B (P = 0.61; U Mann- 
Whitney's test) nor for STAT5A (P = 0.81; U Mann- Whitney's 
test). 

STAT5B immunoexpression levels were significandy different 
between samples grouped according to pTNM classification in 
relation to tumor size (P = 0.045; ANOVA Kruskal-Wallis test), as 
presented in Figure 4. Neuman-Keuls' multiple comparison test 
revealed statistically significant differences between Tla+Tlb and 
T2a+T2b group (see Table 7). 

We didn't find any significant correlations between STAT5B 
immunoexpression level and patients' age (P = 0.61), gender 
(P = 0.40), the history of smoking assessed as PY (P=0.63), and 
AJCC classification (P = 0.17) (ANOVA Kruskal-WaUis test, U 
Mann- Whitney's test followed by Spearman's rank correlation 
coefficient). 
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Figure 4. Box-and-whisker plots, representing STAT5B protein 
levels (ng/ml) in the studied tumor groups classified according 
to pTNM staging in relation to the tumor size. 

doi:1 0.1 371/journal.pone.01 04265.g004 

No significant correlations were found between pSTAT5A 
immunoexpression level and patients' age (P = 0.83), gender 
(P = 0. 1 7), the history of smoking assessed as PY (P = 0.92), tumor 
size according to pTNM classification (P = 0.46), and AJCC 
classification (P = 0.98) (ANOVA Kruskal-Wallis test, U Mann- 
Whitney's test followed by Spearman's rank correlation coeffi- 
cient). 

5. Statistical analysis of the reciprocal relationship 
between STAT5A/B mRNA expression and STAT5A/B 
immunoexpression 

The positive correlation was found between the expression 
levels of STAT5B gene and the immunoexpression levels of 
STAT5B (rho = 0.755, P = 0.04; Spearman's rank correlation 
coefficient) in NSCLC samples, see Figure 5. No significant 
correlation was found between STAT3A gene expression and 
STAT5A protein level (P — 0.43, Spearman's rank correlation 
coefficient). 

Discussion 

It is commonly known that deregulations of STAT signaling 
pathway can result in oncogenesis and the constitutive activation 
of STATS is found in many human malignancies [10,12,24], 
including lung cancer [25]. Our study results confirmed the 
implication of STATS isoforms in NSCLC development. We 
observed the increased expression levels oi STATS genes (A and B 
isoforms) in the majority of studied NSCLC tissues, both at mRNA 
and protein level. Our result are in accordance with observations 
of other investigators who have also reported STAT5A and 
STAT5B overexpression on mRNA and/ or protein level in cancer 
tissue - prostate, breast, colorectal, esophageal and lung - when 
compared to normal tissue [10,24—31]. 

Interestingly, it is known that STAT5A and STAT5B isoforms 
have distinct biological functions [32-33] and the expression 
patterns of these proteins differ among many types of cancers cells. 
Both STAT5A and STAT5B are key tumor growth stimulators, 
but STAT5B overexpression seems to be more significant in head 
and neck and prostate cancers, while STAT5A — in breast cancer 
[31,34,35]. However, the exact role of STAT5A as opposed to 
STAT5B in mahgnant cell transformation is not well understood. 



6 



August 2014 | Volume 9 | Issue 8 | el 04265 



STATS, COX-2 and PIAS3 Expression in Lung Cancer 



Table 7. The results of Neuman-Keuls' multiple comparison test regarding differences in STAT5B immunoexpression levels in the 
studied tumor groups classified according to pTNM staging in relation to the tumor size. 



T1a+T1b T2a+T2b T3+T4 

TIa+TIb 0.0046134 0.92834 

T2a+T2b 0.0046134 0.107350 

T3+T4 0.92834 0.107350 



doi:l 0.1 371 /journal.pone.Ol 04265.t007 

Recently, Tang et al. [36] have demon.strated difTerential 
regulation of mammary carcinoma cell behavior upon forced 
expression of STAT5A and STAT5B: STAT5B exhibited lower 
potency than STAT5A in enhancing survival and anchorage- 
independent growth of mammary carcinoma cells and exerted no 
effects on cell motility. In our study we assessed the immunoex- 
pression of STAT5A and STAT5B proteins in lung cancer tissue 
and found it significantly higher than in normal lung tissue. The 
study performed by Sanchez-Ceja et al. [25] indicated the highest 
nuclear STATS expression in LCC (large cell carcinoma), lower in 
sec and the lowest in AC (adenocarcinoma). In our study we 
included LCC and AC in one NSCC group due to the small 
number of LCC samples. However, we didn't find any significant 
differences between the groups. STAT5A and STAT5B protein 
levels were similar in the studied NSCLC subtypes, however, for 
procedural reasons, we could not analyze the location of proteins 
(nuclear or cytoplasmic). 

The analysis of both gene isoforms on mRNA level in NSCLC 
samples revealed significantly higher expression of STAT5B 
isoform as compared with STAT5A. It might suppose different 
functions for both isoforms also in lung cancer cells, i.e., in result, 
different downstream target genes (including oncogenes) specifi- 
cally induced by either STAT5A or STAT5B. The mechanism of 
lower STATS A gene expression in NSCLC is not known, however 
in lymphomas STAT 5 A was found to be epigenetically silenced 
[37] . Consistent with the obtained results indicating high STAT5B 
expression in lung tumors was our analysis regarding STAT5B 
immunoexpression, as we observed positive correlation between 
these two parameters. 





















4P 




















O CC03 
































































































O 0 










h 

s ° 






O 





0123456789 10 11 
RQ 

Figure 5. Positive correlation between STAT5B protein levels 
(ng/ml) and gene expression levels (RQ values) in NSCLC 
samples. 

doi:1 0.1 371 /journal.pone.Ol 04265.g005 



We documented also statistically significant differences in 
STAT5B expression levels depending on tumors size, i.e., 
significandy higher STAT5B expression in T2 group, which was 
additionally reflected by significandy higher STAT5B protein 
levels in this group. Thus, it could further suggest that STAT5B 
has more important implications than STAT5A in the pathogen- 
esis of NSCLC. As so far there is no available published data 
concerning association of STAT5B expression with tumor 
progression and prognosis in NSCLC. The analysis of STATS 
immunohistochemical staining in NSCLC tissue, performed by 
Sanchez-Ceja et al. [25], suggested the association of STAT5 
protein with advanced lung cancer stages. Regarding other tumor 
types the data are conflicting. The significant association of 
STAT5B protein expression with TNM stage was found in 
colorectal cancer [38], but not in gastric cancer [39]. Interesting 
results were obtained by Tang et al. [36], who proved a dual role 
of STAT5A in human mammary carcinoma cells: promotion of 
tumor formation by enhancing survival and anchorage-indepen- 
dent growth but simultaneous inhibition of tumor metastasis by 
suppression of cellular invasiveness. 

In our study, we tried to estabfish - to our knowledge as a first 
research group - the relationship between STAT5A and STAT5B 
niRNA expression and patient's age, gender, and NSCLC 
histopathological subtype. Such an association could indicate an 
importance o{ STATS mRNA expression level as a diagnostic tool 
in NSCLC. However, we didn't find any significant associations. 

PIAS3 (protein inhibitor of activated STATS), was originally 
identified as a specific inhibitor of the STAT3 signaling pathway 
[40]. Now it is well recognized that the PIAS family (PIASl, 
PIAS3, PIASx and PIASy) regulates a variety of cytokine and 
transcription factors with downstream alterations in apoptosis, 
angiogenesis and a number of signaling pathways. Regarding 
PIAS3 and its role in lung tumorigenesis, it has been demonstrated 
that PLAS3 decreases lung cancer growth and increases the 
antitumor effects of EGFR inhibitors [20] , activates the intrinsic 
apoptosis pathway via altered expression of Bcl-2 and Akt family 
members [21,23], and interacts with several other transcription 
factors, including: ETS, EGR, NR1I2 and GATAl [21], all of 
them playing important roles in cancer development. The STAT3 
independent effect of PIAS3 includes also the repression of 
STAT5 transcriptional activity, as demonstrated by Rycyzyn and 
Clevenger [19]. In our study the expression level ot PIAS3 was 
negatively correlated with the expression level of STAT5B in the 
whole group of studied tumors, and especially, the negative 
correlation was observed in T2 tumor group. 

In several cancers a decrease or loss of PIAS3 expression was 
demonstrated, indicating it as a protein with putative tumor 
suppressor function [22], but a paucity of studies were focused on 
lung cancer [20,21,23]. The results obtained in NSCLC cell lines 
indicated that overexpression of PIAS3 had growth inhibitory 
effect [20,23]. Our analysis revealed the decreased PIAS3 
expression in the majority of the studied NSCLC samples. 
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Moreover, it was significantly lower in NSCC as compared to 
sec. Our results are similar to those obtained by Kluge et al. [20] 
who found PIAS3 mRNA expression higher in SCC than in AC. 

Cyclooxygenase (COX), is the key enzyme in the biosynthesis of 
the prostanoids and plays a central role in many important cellular 
processes, including inflammatory response, tumorigenesis, and 
tumor progression [41]. COX is comprised of three categories, 
including COX-1, COX-2 and COX-3 [15]. COX-2 converts 
arachidonic acid to bioactive lipids including prostaglandin E2 
(PGE2), and their role in initiating and progressing tumors in situ 
was established. It was suggested that the overexpression of COX- 
2 and the resulting increase in PGE2 levels may represent a tumor 
strategy to escape immunosurveillance [42]. CC)X-2 overexpres- 
sion was found by immunohistochemical method in well and 
moderately differentiated carcinomas of the lung, colon, and 
breast [18,43]. In lung cancer, COX-2 overexpression was 
reported to inhibit apoptosis [16], promote angiogenesis [44] 
and metastasis [4.5]. In our study we confirmed intTeased 
expression of COX-2 gene in the majority of NSCLC samples, 
regardless the histopathological subtype. In contrast, it was found 
that COX-2 protein expression in adenocarcinoma was signifi- 
candy higher than that in squamous cell carcinoma [43,46,47]. 
These differences might be caused by further posttranslational 
modifications of COX-2 protein in lung cancer cells. 

Interestingly, we observed statistically significant difiFerences in 
COX-2 expression with respect to tumor size - significantly higher 
gene expression levels were shown in T2a-I-T2b tumors as 
compared with Tla+Tlb group. The meta-analysis combining 
14 published studies, revealed a significant association between 
COX-2 expression and poor survival for stage I NSCLC (i.e., 
encompassing Tla, Tib and T2a tumors) [48]. However, the 
studies reporting the relationship between COX-2 expression and 
survival among lung cancer patients are inconsistent. The study 
performed by Groen et al. [49] confirmed that high COX-2 
expression predicted poor prognosis in NSCLC, but the results of 
more recent study have suggested that the prognostic role of 
COX-2 in NSCLC needs to be confirmed by further high-quality 
prospective studies [50]. 

In human lung adenocarcinoma cells it was found that COX-2 
expression was stimulated via the activation of the STAT5 
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